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An atom interferometer based on a Stern-Gerlach beam splitter is proposed. Atom scattering 
from a combination of magnetic quadrupole and homogeneous magnetic fields is considered. Using 
Raman transitions, atoms are coherently excited into and de-excited from sublevels having nonzero 
magnetic quantum numbers. The spatial regions in which the atoms are in such sublevels are small 
and have magnetic fields designed to have constant gradients. Therefore, the atoms experience 
position-independent accelerations, and the aberration of the coherently separated and recombined 
atomic beams remains small. We find that because of these properties it is possible to envision an 
apparatus producing atomic density gratings with nm-scale periods and large contrasts over 10 — 100 
/im. We use a new method of describing the atomic interaction with a pulsed spatially homogeneous 
field. In our detailed analysis, we calculate corrections caused by the non-linear part of the potential 
and the finite value of the de-Broglie wave length. The chromatic aberration and the effects of an 
angular beam divergence are analyzed, and optimal conditions for an experimental demonstration 
of the technique are identified. 
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I. INTRODUCTION 

An important application of atom interference Q] is the 
production of a periodic spatial profile of the atomic den- 
sity. When an atomic beam propagating along the x axis 
passes through a system of counterpropagating resonant 
optical fields having wavelength A or through a microfab- 
ricated structure having period A/2 directed along the z 
(transverse) axis, an initial atomic state having trans- 
verse momentum p splits into a series of states having 
momenta p + rihk, where n are integers and k — 2tt/X. 
Interference between these momentum states results in 
an atomic density pattern, referred to as a grating, hav- 
ing a period A/2. A detailed review and bibliography of 
the different regimes for producing atomic gratings can 
be found in our recent article B. 

For the purpose of this paper, it is important to un- 
derline that gratings with a sinusoidal density profile and 
nanometer scale period X g <C A, say 



An 



10- 100 ram, 



(1) 



are of particular interest. One method to achieve this 
goal is to use a large angle beam splitter (LABS), which 
splits the initial atomic state into two states having mo- 
menta p ± Ap, where 



Ap = irh/Xg 



(2) 



It was expected that triangular potentials, which one 
can produce with some accuracy using a strong standing 
wave field p|, a magneto-optical scheme [ra, or bichro- 
matic fields Q, can work as a LABS. However, we re- 
cently showed that, asymptotically, the atom den- 
sity profile scattered from this LABS becomes a spatially 



inhomogeneous sinusoid with a period of order X g super- 
imposed on sharp density peaks, separated by A/2. The 
undesired inhomogeneity results from the splitting of the 
atomic state into two groups of momentum states rather 
than into two well defined momentum states. 

In this article we propose a LABS based on the Stern- 
Gerlach beam splitter Our main objective is to 
test the idea of using this LABS to create sinusoidal 
atomic gratings with periods much smaller than the op- 
tical wavelength but with coherence lengths much larger 
than the optical wavelength. 

The article is arranged as follows. In the next Section, 
we outline the principles of operation of the proposed 
devices. In our detailed calculations, we then consider 
atomic scattering from a finite thickness layer of a linear 
potential (Section III). Corrections for weak and strong 
acceleration are evaluated in Sections IV and V. The ap- 
pendix |A| is devoted to calculations of the Raman tran- 
sitions between Zeeman sublevels. We summarize our 
results in Sec. VI, and discuss their applicability to a 
beam of 87 Rb atoms. 



II. PRINCIPLE OF OPERATION 



From a naive point of view it seems impossible to 
use Stern-Gerlach beam splitters for atom interferome- 
try, since such beam splitters produce inerferometer arms 
corresponding to different Zeeman sublevels. These arms 
could not interfere. This problem can be avoided as fol- 
lows (see Fig. [l]) 
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FIG. 1. Principal scheme of the atom interferometer based 
on the Stern-Gerlach beam-splitter. 

In Fig. ^ we show the principle of the atom interfer- 
ometer of interest. If an atom having angular moment 
G = 1/2 is polarized along the x or y axis, then after 
splitting by the Stern-Gerlach magnet I, the two arms of 
the interferometer contain atoms in orthogonal but mu- 
tually coherent states having Zeeman quantum numbers 
rn = ±1/2, respectively. Magnets II and III reflect the 
atomic beams in the two interferometer arms, and re- 
combine them. Before the recombination, the internal 
state in one arm of the interferometer is flipped by ap- 
plying an RF-induced or optically induced 7r-pulse in a 
spatially homogeneous magnetic field. All atoms arrive 
in the recombination region in the same internal state, 
and interference occurs. 

An ideal sinusoidal grating would arise if an incident 
plane atomic wavefunction is split into two coherent com- 
ponents in a spatially homogeneous magnetic field gradi- 
ent. Evidently, a realistic Stern-Gerlach magnet, charac- 
terized by a number of higher order spatial derivatives of 
the magnetic field, does not satisfy this requirement. To 
realize a magnetic field with constant field gradient, we 
consider the superposition of a homogeneous bias field 
B s directed along z axis and a quadrupole magnetic field 
produced by four (anti-)parallel currents propagating in 
the ±?/-directions. This combination allows us to real- 
ize a small volume around the center of the quadrupole 
field in which the field gradient is approximately con- 
stant, with no zero of the magnetic field being present. 
Zeroes of the magnetic field need to be avoided in or- 
der to prevent non-adiabatic spin flips (Majorana tran- 
sitions). Since there are no zeroes of B, an atom in a 
magnetic sublevel m with respect to a quantization axis 
identical to the B-field direction adiabatically follows di- 
rection changes of the _B-field along the atom's trajectory, 
and will always remain in that sublevel. The bias field 
also allowed us to generate potential with a dominant 
linear term and small nonlinear corrections. The poten- 
tial for the atom's center-of-mass motion is then given 
by V(r) — gfiBfn |B(r)|. Atoms prepared in the m = 



sublevel do not interact with the field at all. 

Based on the preceding considerations, we can now 
outline how we realize Stern-Gerlach beam splitters with 
small aberration. We assume that the above described 
combination of £>-fields produces a region of almost con- 
stant field gradient centered at the origin. A plane matter 
wave in state m = propagates in the -t-x-direction. The 
wave is not refracted or diffracted by magnetic-dipole 
forces while it traverses the fringe fields of the magnets. 
External coupling fields are applied in a narrow plane 
located at x — — d, splitting the wave into a coherent su- 
perposition of two waves in different (relevant) Zeeman 
sublevels (e. g. m — ±1),. The spin components of the 
wave experience an acceleration oc to, and acquire a dif- 
ferential transverse momentum change Ap. A second set 
of coupling fields in a plane at x = +d return the ac- 
celerated atoms into the state m = and terminate the 
interaction with the magnetic field. As a result, the sys- 
tem coherently splits an atomic plane wave entering in a 
single magnetic sublevel into two momentum components 
exiting in the same magnetic sublevel. aberration effects, 
i.e. unwanted curvatures in the phase fronts of the out- 
going waves, are minimized by the fact that all magnetic 
acceleration is localized to a small region — d < x < d, in 
which the field is not strongly contaminated with higher- 
order multipole terms ("fringe fields"). 

If the momentum states overlap spatially, an atomic 
grating will form. The splitting between Zeeman sub- 
levels caused by the external coupling fields determines 
the momentum space distribution and the properties of 
the atom grating in the detection plane. The grating 
phase is determined by the difference between the atomic 
wave functions phases acquired along the two arms of 
the interferometer. Owing to the phase sensitivity to 
the atom velocity and the magnetic field instability, the 
grating can be washed out, or one has to require that the 
atomic beam and field characteristics must be beyond 
the current state-of-art. One has to choose a splitting 
scheme and interferometer geometry that minimizes this 
sensitivity. Two types of interferometers are shown in 
Fig. | 

a b 
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FIG. 2. Schemes to create atomic gratings with (a) asym- 
metric and (b) symmetric atom interferometers. I, II, and III 
are beam-splitters, whose possible layout is shown in Fig. |^. 

In the triangular interferometer (Fig. Ha) atoms in 
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arms 1 and 2 have different kinetic energy The phase, as- 
sociated with this difference, the so called Talbot phase, 
leads to periodic oscillations of the density distribution 
||, initially discovered for light |l0) and also observed 
in an atom interferometer pl| . The Talbot phase plays 
a critical role in time-domain atom interferometry Jl2"| , 
where it has been used for precise recoil frequency mea- 
surements pl|. For interferometers in the spatial do- 
main, the Talbot phase degrades the atom grating, and 
one should prefer a symmetric interferometer (Fig. ^|b) 
that produces no phase difference. 

To produce a symmetric interferometer one first has 
to split an initial atomic state \m = 0) symmetrically 
between, for example, states \m = ±1) . Co-propagating 
cross-polarized optical waves can produce two-quantum 
transitions between Zeeman sublevels via excited state 
manifold to create ground state interferometers. Our cal- 
culations show that a magnetic field in the acceleration 
zone can lead to a Zeeman splitting u>z larger than the 
inverse interaction time r^, 

uzn > l, (3) 

and effective coupling occurs only if one superposes two 
waves at frequencies fi and SI + ujz needed for a two- 
quantum resonances. Next, the Zeeman splitting is typ- 
ically comparable with the excited state hyperfine split- 
ting uJhf, so that, rigorously speaking, to obtain two- 
quantum transition amplitudes one has to know the ex- 
cited state manifold structure for a given magnetic field. 
The situation simplifies if the detuning A between the 
waves' frequencies and the ground-excited state transi- 
tion frequency is larger than both the Zeeman and hy- 
perfine splittings, i.e. 

|A| > max{v z ,0Jhf}- (4) 

In this case the reduced matrix element of the two- 
quantum transition is proportional to that arising in the 
absence of Zeeman and hyperfine splittings. Since alkali 
ground states have angular moment Jq = 1/2, selection 
rules allow field absorption and emission processes only 
where the angular moment projection changes by at most 
one, and thus the two-quantum reduced matrix elements 
between ground state levels vanish. As a result, a combi- 
nation of <r_|_ — <r_ fields produce no transitions. For this 
reason we assume here that the fields' frequencies and 
polarization vectors are chosen as {Jl,z}, {f2 + uiz, x}, 
while the atomic initial state is \G — 1, m = 0) , where G 
is the total atomic angular moment (see Fig. IJd). In the 
atomic rest frame traveling waves localized in the thin 
layer act as a pulse. For G = 1, a proper choice of the 
pulse area allows one to split all atoms symmetrically be- 
tween Zeeman sublevels |m = ±1) (see appendix |A|) and 
start their acceleration. 




m = - 1 

FIG. 3. (a) Scheme of the beam-splitter. The beam 
splitter involves four currents Uj comprising a magnetic 
quadrupole, a homogeneous bias magnetic field B s , two 
thin layers of optical fields located at x = id. Initially 
atoms move without acceleration in the \m — 0) state. 
In the first layer, atoms split between Zeeman sublevels 
jm = ±1) and start accelerating. In the second layer atoms 
are partially returned to the initial state and stop ac- 
celerating, while atoms in other sublevels (their trajec- 
tories shown by dashed curves) leave the interferometer. 

(b) Coupling of the atomic Zeeman sublevels by optical waves 
propagating along the y axis inside each layer and having fre- 
quencies and polarization vectors {fi,z} and {fi + ojz,x}. 

To stop the acceleration, one applies another set of 
traveling waves located on a plane x = +d to return 
atoms back to the \m = 0) states. The x-range within 
which the acceleration that split the atomic beam are 
active is thereby limited to the thin region \x\ < d. One 
has to distinguish the regimes of weak and strong accel- 
eration, characterized by 

T) < 1, (5a) 

V > 1, (5b) 

respectively, where 

r) = Az/b, (6) 

b is an atomic beam radius, and Az is the atom displace- 
ment during the acceleration. When the acceleration is 
weak (the case shown in Fig. ^a), on the planes x = —d 
and x = +d the same pair of fields can be used. Since 
the Zeeman splitting is equidistant, these fields drive a 
chain of transitions \m = ±1) — > |m = 0) — * \m = ^fl) , 
and, evidently, can not return all atoms to the \m = 0) 
state. Nevertheless, it is possible to maximize the am- 
plitude of return to the |m = 0) state. At x > d, atoms 
that remain in the state |m = ±1) continue to accelerate 
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and eventually leave the interferometer, while atoms in 
the state \m = 0) are of further interest. 

In Fig. ||, quadrupoles II and III act along the spa- 
tially separated arms of the interferometer, and are ad- 
justed such that they reverse the transverse components 
of the atomic momenta. For this purpose, one may still 
apply fields {f2, z} and {ft + u>z, x} to transfer atoms at 
x = —d to the states |m = ±1) , decelerate them between 
— d < x < d, and return them back to the state \m = 0) 
at x = +d. (the origin is assumed to be at the center 
of the respective quadrupole field). Several additional 
arms of useless atoms in |m = ±1) states will be pro- 
duced. Our calculation show that only one-eighth of the 
atoms will be properly recombined to produce a grating, 
while seven-eighths will be lost. To avoid the loss, we 
propose to use a second hyperfine manifold, having an- 
gular momentum G = 2. If at 2; = — d one applies fields 
{f2, z} and {fl + W2i,i0) x } , where u>2i,io is the frequency 
of the transition \G = 2, m = 1) — > \G = 1, m = 0) , then 
only this two-level scheme is involved, because the fre- 
quencies nuj2i,iQ (n ^ 1 is integer) no longer coincide 
with any atomic transition frequency. Choosing a field 
pulse area A = ir, one can transfer 100% of the atoms 
at x = —d into the \G = 2,m= 1) state, accelerate the 
atoms, and return all of them back to the \G = 1, m = 0) 
state at x — d by another 7r-pulse. 

After the action of quadrupoles I - III of Fig. ||, the 
total difference of momenta in the interferometer arms is 
given by 



Ap = Api + Ap n . 



(7) 



Momentum kicks Api associated with quadrupole i have 
opposite signs. One needs to use Apn partially to cancel 
Api. The larger |Ap//| requires a larger field and more 
severe conditions for the gradient homogeneity. As a re- 
sult, for a given desirable gratings period X g it is better 
to choose 



such that 



|Ap/| < \A PII \ 



Ap H « Ap. 



(8) 



(9) 



In this case the role of quadrupole I is just to split the 
beam into two arms, while quadrupoles II and III are 
responsible almost entirely for the grating formation. 

Owing to the inhomogeneity of the field gradient, the 
finite time of the interaction, and the finite angle of the 
atom scattering, instead of changing the atom momen- 
tum p by fixed value Ap, one produces a wave packet in 
the momentum space near the momentum p + Ap with 
a width that increases for larger momentum kick Ap (or 
smaller \ g ). We analyzed this effect recently for LABS 
produced using resonant fields 0,||. The atom grating 



profile would be damaged if the wave packet width be- 
comes larger than fijb. In this article we evaluate cor- 
rections to the wave function associated with the fac- 
tors listed above. For a given grating period A s , we find 
other important characteristics of the problem from the 
requirement for corrections to be small. These wave func- 
tion corrections allow us to choose an atomic beam aper- 
ture b and velocity u, the length of the interaction zone d, 
the magnetic field gradient B' , and the bias field strength 
B s to obtain a desired grating with a given accuracy. 

The performance of the various Stern-Gerlach accel- 
eration regions in the above schemes is limited by chro- 
matic aberration and aberrations due to inhomogeneities 
of the field gradients. The detailed analysis presented in 
the following sections provides a quantitative foundation 
to estimate these effects, and to identify the best possible 
operating conditions. 



III. ATOMIC SCATTERING FROM A LINEAR 
POTENTIAL WITH SMALL CORRECTIONS 

Our Stern-Gerlach beam splitters can be characterized 
by a potential 



U (z, x) = Uq (x) — f (x) z + Ui (z, x) , 



(10) 



which acts on atoms propagating predominantly in 
the ir-direction. The potential acts in the narrow layer 
\x\ < d, and consists of a homogeneous part Uq (x) , a 
large linear part — / (x) z, and a small nonlinear addition 
Ui (z, x) . In the following, we analyze the propagation of 
matter waves in such a potential. 

A time- independent solution $ (z, x) of the matter 
wave at a given energy E in the potential U (z, x) fol- 
lows the Schrodingcr Equation 

^M+U(z,x) $(z,x)=EQ(z,x), (11) 

where p and M are the momentum operator and the 
atomic mass. When the potential is weak compared to 
the kinetic energy, which is mostly given by the motion 
in x-direction, 



\U (x)\ «£, 



(12) 



one can use a slowly varying amplitude approximation 
for the wave function. Introducing a "time" t = x/u (x) , 
where 

u{x) = [2(E- U (x))/M] 1/2 (13) 
is an atomic velocity, one can seek a solution of the form 



$ (z, x) = exp 



i— ( dxu(x) 
h J-d 



(14) 



where rf> (z, t) is the slowly varying wave function am- 
plitude (referred to below simply as the wave func- 
tion). In the momentum representation, ^ (p, t) 



4 



-1/2 

obeys the equation 

,,2 



(2nh)~ ' J dz exp (—ipz/K) ip (z, t) , this wave function Recently, this equation has been exactly solved in the co- 
ordinate representation Jl^-[l6t . For the purposes of this 
article, we derive the solution with an alternate method 
using the momentum representation and an accelerated 



ihd t V 



P 

2M 



ihf (t)dp + Ui(ihdp,t) 



*, (15) 



frame 



where / (t) — f (ut) is a force, U\ (z, t) = U\ (z, ut) , and 
Q represents small terms arising from a second derivative 
in time and the slow variation of the atomic velocity, 



p = p Q + 5p (t) , 



Sp(t) 



h 



Q=-i-u'(x)[l + i J ^d t )- 



dhf(tx), 



2Mu 2dt 



(16) 



Neglecting U± and Q one arrives at a one dimensional 
Schrodinger equation with a time dependent spatially ho- 
mogeneous force, 

iK(d t + f(t)d p )*(p,t) = ^(p,t) (17) 



where the wave function evolves as 



(18a) 
(18b) 



(19) 



Solving this equation and returning back to the lab frame, 
one finds the following common expression: 



(p, t) = exp 



* dti ^ + 5p{t 1 )-6p(t)] i 



2Mh 



V(p-5p (t) , -r) . 



(20) 



If Eq. (|T^) is written in the accelerated frame (|T^) , the term proportional to p\ is responsible for the matter wave 
spreading. One can neglect this term if Tp\j2MTi <C 1. In the case of a diffraction- limited single-mode atomic beam, 
po is a momentum typical of the atomic beam spread, po ~ Ti/b, where b is the radius of the incident wave. Then, the 
just mentioned condition is equivalent to 



m = T l T s < 1, 



(21) 



where r s = Mb 2 /Ti is a time characteristic of the spreading of the matter wave. Assuming that this condition is valid 
we drop the quadratic term in p , and arrive at the equation 



ihdt^ (p ,t) 
Seeking a solution of the form 



Sp 2 (t) Sp (t) 



2M 



+ 



M 



-Po + Ux (ihd po ,t) + 



*(Po,*) 



where 



* (po, t) = exp [~i(f> - ip Q Sz (t) /h] * (p , t) , 
8p(h) 



6z (t) 



dti 
dti 



M ' 

Sp 2 (h) 
2Mh ' 



one finds that $ (po, t) evolves as 



im = {Ui [Sz (t)+ihd P0) t] + Q}it>, 



where 



Q = -i-v! (x)[l + i 



Mu 2 



2Mu- 



q = exp [i(f> + ip 5z (t) /%] (d t - f (t) d Po ) exp [-i<j> - ip 5z (t) /%] . 
The wave function in coordinate space at the exit of the interaction zone, if) (z, r) , is given by 

dp 



tp (z, t) = exp [~i(f> + izAp/h] 



(2tt?i) 



Yj^ exp [ipo (z - Az) /%] * (p , r) , 



(22) 

(23) 

(24a) 
(24b) 

(25) 

(26a) 
(26b) 

(27) 
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where 



Ap = Sp (r) and Az = Sz (r) (28) 

is a classical change of the atomic momentum and position under the spatially homogeneous acceleration acting for 
a time 2r. One can consider an atomic grating close to sinusoidal, if the period X g is smaller than the transverse 
extension b of the matter wave, which means that 

72 = ft/Ap6<l. (29) 

In the zeroth order approximation in U± and Q, ^o(po,t) — const and, therefore, the zeroth-order wave function 
( p7| ) is given by the expression 

ipo [z, t) = exp [—i<f> + izAp/h] ip(z— Az, — r) , (30) 

which can be obtained also from the common solution ( p0[ ) at the assumption (pl|). One sees that a matter wave moving 
with a sufficiently large and time-independent (to neglect term Q) velocity u through a layer of the homogeneous force 
for a time smaller than the packet spreading time is just displaced in the phase space along the classical trajectory. 

In the absence of the higher-order effects described below, a purely sinusoidal grating can be formed by interfering 
atomic momentum components. 



IV. WEAK ACCELERATION 



In the following two sections we calculate higher order effects that degrade the ideal scattering behavior of the 
matter wave. 

Consider an atomic beam propagating with velocity u along the x-axis and interacting with a quadrupole magnetic 
field B (z,x) produced by four currents, directed along the y-axis, located in the (z,x) plane at (ia,ja x ) = ±1), 
and given by Jjj — lij. In addition to the quadrupole field, one applies a spatially homogeneous magnetic bias field 
B s = zB s , such that the total magnetic field is given by 



B z (z,x) 
B x (z,x) 



-B }b s +a ^2 i3{x-ja x )j (z - ia) 2 + (x - ja x ) 2 > 



Br, a 



ij (z — ia) / (z — ia) 2 + (x — ja x ) 2 



(31a) 
(31b) 



where Br> is the absolute value of the magnetic field of one current at a distance a and b s — —B s /Bq. In the rest 
frame (x = ut) , an atom in the internal state characterized by orbital angular momentum Lq, electronic spin total 
electronic angular momentum Jq, nuclear spin /, total angular momentum G, and projection m of the total angular 
momentum on the magnetic field direction, moves in a potential 



where 



fj, = {1 + [J G ( J G ■ 
x {[G(G + 1) 



U(z,t) = fj,\B(z,ut)\, 



1)-L(L + 1) + S(S+1)}/ [2J G (J G + 1)]} 
I- Jg (Jg + 1) - I (I + 1)] / [2G (G + 1)]} (i B m 



(32) 



(33) 



is the projection of the total magnetic moment on the direction of B, and fig the Bohr magneton. 

We assume that the atomic beam is centered at z — and has a radius b -C min{a, a x } , and that the ^- or ir- 
Raman fields, which turn the interaction with the magnetic field on and off, are located at x — ±d (d min{a, a x }) , 
such that the half-duration of the interaction with the potential ( |32"| ) is r = d/u. When b s ^ 0, one can expand the 
potential ( |32"| ) in the vicinity of the (z = 0,t = 0) point. Omitting homogeneous term of Eq. [l(] and assuming, for 
simplicity, that there is no explicit time dependence of the force, one finds 



U(z,t) = -fz+Ui [z,t) 



OO OO 



EE 

n—2 rn—0 



Ui (z,t) = f 2^ 2^ Cnm 



(a,b s ) 



(34a) 
(34b) 



G 



where 

/ = (35) 
and B' is the magnetic field gradient at the quadrupole center, and r a — a/u 7 . The dimensionless coefficients 

c nm (a, b s ) = /n , (2 ^ }! W2m ^ t=Q (36) 

depend on the quadrupole size ratio 

a = a x /a (37) 



and the bias field's relative strength b s . For a magnetic quadrupole, the coefficients c nm , used in further calculations, 
are given by 



C21 = 16a 
c 2 2 = -16a 



16a 2 - 36 2 (a 2 - l) (a 2 + l) 2 j b; 3 (a 2 + l)~ 6 ; 



(38a) 

1536a 4 - 352a 2 6 2 (a 2 - l) (a 2 + if + 6 4 (a 2 + l) 4 (21 - 62a 2 + 21a 4 )] 6J 5 (a 2 + l)~ 10 , (38b) 
16a |-122880a 6 + 353285 2 a 4 (a 2 - l) (a 2 + l) 2 - 326 4 a 2 (a 2 + l) 4 (95 - 234a 2 + 95a 4 ) 
+ b G s (a 2 + l) 6 [81 (a 6 - 1) - 463a 2 (a 2 - l)] } b~ 7 (a 2 + l)~ 14 , (38c) 
c 30 = 2(a 2 -l) (a 2 + l)- 2 . 



f'23 



To indicate the structure of the c nm , in the following 
matrix indices (n, m) with c nm ^ are marked by a "V" : 



(Aid) 
(43) 





TO = 


1 


2 


3 


n = 2 




V 


V 


V 


3 


V 


V 


V 


V 


4 




V 


V 


V 


5 


V 


V 


V 


V 



(39) 



rj = K89/2/3, 

while the small parameters (Efi, E9f) are given by 



m = (27r)- 1 69e g (3- 2 , (44a) 



(44b) 



We characterize the problem by dimensionless param- 
eters 



= b/a, 


(40a) 


S = d/a, 


(40b) 


Eg = Xg/a, 


(40c) 


= XdB/Xg, 


(40d) 


X dB = 2irh/Mu 


(41) 



One can express the atomic beam and magnetic 
quadrupole characteristics through the parameters (fill). 
For example, using Eqs. ( gog , £0]) and then @ p5]pb[), 
one finds the atom velocity and magnetic field gradient: 



u = (2Trh/MXg)0- 1 , 
B' = (mr 2 h 2 /nMXf) e^ -1 * -1 . 



(45a) 
(45b) 



where 



is the atomic de-Broglie wavelength and 9 is the angle of 
atom scattering. The constant force one needs to apply 
to achieve a given atomic grating period can be found 
from Eqs. (§| @ |) to be 



One can use Eq. (pq ) to calculate corrections to the 
unperturbed atom wave function tp Q (po,t) = -0 (po, ~ T ) ■ 
Using the estimate 



hi:), 



/'o 



(46) 



/ = mtTiu / 2dX Q 



(42) 



where the parameter k sa 1 for quadrupoles II and 
and k = Api / Ap <C 1 for quadrupole / (see Fig. ||b). 
Consequently, the atom displacement Az = ndO/2 and 
the parameter (ph is given by 



for the case of weak acceleration (f3a|), one can neglect 
the term 5z (t) in Eq. (25). After this, one can calculate 
the first-order correction ^i( pq,t ), associated with the 
(n, to) term of the expansion (|24b| ) . Substituting the ex- 
pression for ^i(pq, t) in Eq. (|27|) one finds the correction 
in the coordinate representation 
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fa (z, t) = -i [2fc nm d 2m+l I (2m + 1) ha n + 2m - x u\ (z - Az) n (z, r) , (47) 

where ipo (z, r) is given by Eq. ( |3C| ) 

We now proceed to calculate the correction ^ q (po,T) associated with a term Q in Eq. (25). In Appendix [b] we 

found the conditions under which one can neglect the first term in the Eq. (26a), while the operator (26b) is reduced 
to the expression 

q~-fd po . (48) 

Consequently, the relevant correction in coordinate space is given by 

^ q (z, r) = -i (fr/fiMu 2 ) {z - Az) 2 % {z, t) . (49) 

In contrast to (f47|), this correction arises from the wave packet motion through a field with a homogeneous gradient. 

The parameters f3, S and 9 have to be chosen such that corrections (^7], fl9]) are small. We estimate these corrections 
at z — Az — b. Introducing a small parameter 

e q = \% (Az + b, t) /V> (Az + b,r)\ (50) 

one finds for 9 

9 = e q e g e q S(3- 2 , (51) 

where 

e q =8n- 1 K - 2 (52) 

Among the corrections tp± for different n and to, the leading terms arise from 

(n,m) =(2, mi), (53a) 
(n,m) = (n 2 ,0) , (53b) 

as they are the first non-zero terms in the row n = 2 or column m = of the matrix (|39|), and the other terms are 
higher powers of the small parameters r/r a = 6 or \ (z — Az) /a\ % p. Introducing the corresponding small parameters 
£i and e 2 for the relative weight of corrections related to Eqs. ( |53| ) and using Eq. (^), one arrives at equations 

§ 2m lj3 2 _ eiEgSu ( 54a ) 

/3" 2 = e 2 e g e 2 , (54b) 
ei = (2toi + 1) (kit \c 2mi ly 1 , e 2 = (ktt |c„ 2 o|) _1 ■ (55) 



where 



Solving Eqs. (|54| ) one finds 

s = £ y 2 ™ l£ - i /™i«2 £ (n 2 -2)/ 2mi „ 2/(5 (a)6s) ; (56a) 

f s (a,b a ) = e ; /2mi e 2 - 1/mi " 2 , (56b) 

/^(e^) 17 " 2 /^",^), (56c) 

/^(a,6.) -e 2 /n2 (56d) 

= £ y 2 ™i £2 - (2mi+1)/mi " 2 £ g 4 + (" 2 - 4mi - 2 )/ 2mi " 2 /9 (a, &,) , (56e) 

fe (a, b s ) = e\ /2mi e 2 (2mi+1)/min2 e q . (56f) 
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One can use Eqs. (|56j) to estimate the efficiency of the 
beam splitter with arbitrary non-linearity. 

We now return to the quadrupole configuration at 
hand. The non-linearities depend on the ratio of the 
quadrupole sizes a and the relative strength of the bias 
magnetic field b s . These quantities can be used to dimin- 
ish the role of the non-linearities. One can choose them 
such that either the coefficient C21 or the coefficient C30 
vanishes. Our calculations show that it is more effective 
to choose C21 = 0; therefore, in the remainder of this 
article, we consider that case. Equation 

021 = (57) 

determines the ratio of the size a as a function of b s . The 
function a (b s )\ c _ is shown in Fig. ^a. 



linearity that is quadratic in space and 4th-order in time 
also vanishes, i. e. 

c 22 = 0. (58) 

So, one has to consider the next term, (2,3), in the 
row n = 2 of the Table ( |39|). When mi = ri 2 = 3 and 
(a, b s ) is a root of Eqs. (|57|, j58| ), one finds from Eqs. 

6 = 1.45 £ J /6 e2 1/9 4 /18 5 ( 59a ) 
(3 = 0.871 (e 2 e 3 ) 1/3 , (59b) 
9 = 3.19Kr 25 / 18 £l /6 £ - 7/9 e9 4/ 18 . (59c) 
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FIG. 4. Dependence of the quadrupole size ratio (a) , 
width of the acceleration zone (b) , atomic beam radius (c) , 
and scattering angle (d) on the relative strength b s of bias 
magnetic field. 

For c 2 i = the lowest-order correction is the term 
corresponding to mi = 2, n 2 = 3 (see matrix (|39|)). The 
functions f[a(b s ) ,b s ] are plotted in Fig. [| for k = 1. 
These functions have notable features for the relative field 
strength & s « 0.537, where the ratio of quadrupole size 
([37]) is a w 1.99. These singular features arise because 
at & s w 0.537 - in addition to condition fl57D - the non- 



V. STRONG ACCELERATION. 



When rj ^> 1, one can expand the nonlinear part of 
the potential U\ in Eq. ( |25| ) in the operator ihd Po . The 
zeroth-order term Ui [Sz (t) , t] depends only on time and, 
therefore, changes the phase of the atomic wave function 
(p3|) to the value 



5p 2 (t0 



2M 



+ Ui[8z (t0,*i 



(60) 



After a phase transformation one arrives at the equation 



im = {u z {t) ihd P0 + Q}v, 



where 



u z {t) 



dUi{z,t) 



dz 



(61) 



(62) 



z=Sz(t) 



Since one can still neglect the first term in Eq. I26a|) 
and use Eq. (|i^) for the operator q (see Appendix ||), 
the previously calculated correction tp Q (z, r) , associated 
with the Q - term, and, therefore, Eq. (|5l| ) are still valid. 
Using the expansion ( 34b ) one obtains a series for the op- 
erator ([32j). Keeping only the (n, m) term of this series 
in the right-hand-side of Eq. ( |6l| ) , one finds that the cor- 
responding correction in the coordinate representation is 
given by 



l\)l (z,t) 



nc nm I nm pd 2{ - m+ ^- x {z - Az) /ha 



2m+n- 



l u 2n - 



- 1 (2M)™- 1 1 Vo (z,t) 



where 



Inm = y (1 + 



2(1-1) 



(63) 



(64) 



Requiring the magnitude of this correction at z = Az + b to be £1- and £ 2 -times smaller than the zeroth-order solution 
■00 (z,t) for (n, m) = (2, mi) and (n 2 ,0) , respectively, and expressing the atom velocity u and the force / through 
parameters (40), one obtains equations 
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Sl+2mi q 

\T12 — 1 



PS 1 -" L v = ei£g£x, 



13(59) " " = e 2 e g e 2 , 
where for this case we define parameters 6j as 

ei = 8/K 2 nC2 mi h mi , e 2 = 2 3 " 2 ~ 2 / K " 27m 2C, l2 o^ 2 o- 
Equations ( |5l| , |65|) constitute a system for three variables S, (3, 8, which has the solution 

S = e pn 2 -3)/ 2 , e -X/2 7 (eg/ ^ )( „ 2 - 2 )/ 27 ^ ( ^ 5 



(2n 2 -3)/2 7 -l/2 7 ( n2 _2)/2 7 
e l e 2 e g > 

_ ("2-l)/ 7 -(mi+l)/ 7 (n 2 -2)(mi+l)/ 7 mi(n 2 -l)/ 7 



/3 = 4' 



//3 (a,6 s ) 



fe (a, & s 



_ (2n 2 -l)/2 7 (4mi+3)/2 7 / , n [2(m 1 + l)-n 2 ]/2 7 , , , s 
— £l £ 2 /el/Ms]; 

_ _(2n 2 -l)/2 7 (4mi+3)/2 7 -[2(m 1 + l)-n 2 ]/2 7 



/3 = 1.47«- 2 / 5 el /5 (e g /e 2 )^ 
= 1.70049k- 5/4 £2 /4 (e s /£i£ 9 ) 1/4 



where e g is defined by Eq. (|52j) and 7 = n 2 — 2 + mi (2n 2 — 3). 

The further consideration is the same as in the previous 
Section. If one chooses the quadrupole axes ratio a such 
that the potential's quadratic term in space and time 
vanishes, i.e. if a = a (b s ) , where the function a (b s ) is 
defined explicitly by Eq. (|^) and shown in Fig. ^a, then 
again mi = 2 and n 2 = 3. The functions / [a (b s ) , b s ] are 
shown in Fig. H. 



^2/53/5 

£ q 1 



(65a) 
(65b) 



(66) 

(67a) 
(67b) 
(67c) 
(67d) 
(67e) 
(67f) 

(68b) 
(68c) 



VI. DISCUSSION 




FIG. 5. Strong acceleration regime. Dependences of the 
width of the acceleration zone (a) , atomic beam radius (6) , 
and scattering angle (c) on the bias magnetic field relative 
strength b s . 

At the point of divergence, b s ss 0.537, choosing mi = 
n 2 = 3, one finds 



8 = lAAn-^e\' 2 \e a /e 2 e q )^\ 



(68a) 



The scattering of an atomic center-of-mass motion 
wave packet from a narrow layer of quadrupole and 
bias magnetic fields is analyzed. The combination of 
these fields produces an approximately linear potential 
for atoms. It was shown that for a purely linear poten- 
tial, infinitely small atomic de-Broglie wave length, and 
time of interaction smaller than the wave packet spread- 
ing time, the wave packet scatters along the classical tra- 
jectory, changing its momentum by a given amount Ap 
without any wave packet deformation. Only when this 
regime of scattering is realized, at least approximately, 
one can expect that an interference between scattered 
and recombined components of the atomic wave function 
leads to a sinusoidal atomic grating of nanometer-scale 
periodicity. 

In this paper we calculated corrections to the atomic 
wave function caused by potential non-linearities and a 
small atomic de-Broglie wave length. When the grating 
period X g , the quadrupole size a along the grating forma- 
tion direction, the relative strength of the bias magnetic 
field b s , and the relative weight of the corrections caused 
by non-linearities (e\ and e 2 ) and by a finite de-Broglie 
wave length (e q ) are given, one can use our analysis to de- 
termine the atomic beam velocity u and transverse size 6, 
the magnetic field gradient B'&nd the bias field strength 
B s , the thickness of the interaction layer d, and the ratio 
of the quadrupole axes a that will minimize nonlinear 
effects. 
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One can consider nonlinear corrections to the wave 
function as a spherical aberration of the beam splitter. 
There are two more types of aberrations, namely chro- 
matic aberration and the atomic beam angular diver- 
gence. Chromatic aberration arises from averaging the 
grating over the atomic longitudinal velocity u. Since the 
atom momentum change Ap is proportional to the time 
of acceleration r oc l/u, the grating period [by Eq. (|^)] 
is linear in u, 



X g oc u. 



(69) 



To achieve a grating with a given period, one has to use 
a monovelocity beam. Chromatic aberration occurs as a 
consequence of a small but finite width of the velocity 
distribution. To be specific, consider a Gaussian distri- 
bution 



W(u) 



r-V2 



exp 



(70) 



where uq is a mean velocity, a — (u — uq) /uq, and, 
for a beam having longitudinal temperature O, Co = 
u^ 1 (2fcs9/M) 1 ^ 2 is the small relative width of the dis- 
tribution. For a given velocity, atom interference results 
in a term in the atom density, 



p oc cos (4>tot) 



where iptot is a total phase difference of the wave func- 
tions in two arms of an interferometer. The sensitivity 
of the grating to the velocity results from the velocity 
dependence of ifitot- Our purpose was to create phase 
difference 



k g z : 



1 



(72) 



where <p z o is the phase at u = uq, k g = 2ir /X g is a wave 
number associated with the grating period, and we take 
into account Eq. (|69|). 

The largest phase that the atoms acquire is the Talbot 
phase associated with the atomic kinetic energy. This 
phase leads to Talbot oscillations |)) of the interference 
pattern, first observed in a Na beam . Averaging over 



the longitudinal velocity is equivalent to the averaging 
over the Talbot phase. One has to choose an interfero- 
metric scheme, in which Talbot phases can be compen- 
sated. Evidently, the symmetric configuration of the in- 
terferometer satisfies this requirement. Moreover, for this 
configuration one compensates not only Talbot phases 
acquired during the free particles propagation but also 
those associated with the atoms' acceleration inside the 
beam splitters, independently of the acceleration time. 

The next contribution to 4>tot, which we denote as (j)u, 
is caused by the fact that accelerated atom has slightly 
different velocity u (x) during acceleration, owing to the 
homogeneous part of the potential Uq (x) [see phase fac- 
tor in Eq. ([[J)]. Since, during acceleration in the sym- 
metric configuration, atoms in two arms are in substates 
having opposite magnetic quantum numbers, they ac- 
quire phases of the same magnitude and opposite sign. 
Therefore, the grating phase is twice as large as the phase 
along a given arm. For the arm 1, expanding Eq. (Jl^) to 
first order in Uq (x) / E, one finds 



-d 



dx 



U (x) 



(73) 



~ 4>ua (l - o + a 2 ) . 
(71) For the potential produced by the magnetic quadrupole, 



This phase behaves as l/u, i.e. 



U (x) « fab s (l + a 2 ) /8a, and therefore 
0(70 = Kb s (l + a 2 ) 2 /4as g . 



(74) 



In the case of weak acceleration, there are no other con- 
tributions to the phase associated with chromatic aber- 
ration. Owing to the large value of this phase, even for 
small widths of the velocity distribution, the grating can 
be washed out after averaging over velocities. To avoid 
this situation, we propose to insert one more element into 
the interferometer, a region of homogeneous magnetic 
field. In this region an atomic wave function acquires 
a phase <j) c oc u , i.e., 4> c = 4> c q (l — a + of) . Choosing 
this phase to compensate the phase <fiu (|73j), one finds 
for the grating averaged over velocities 



-1/4 



{kza f /4 (1 + C 2 ) cos|/c 



;z — 0.5 tan 



^(C) + (fc^o) 2 C/4(i + C 2 )}, 



(75) 



p oc (1 + C ) exp 
where £ = kza 2 . When 

C « 1, (76) 
we see that averaging over velocity creates a Gaussian envelope of the grating profile having a half width 

s = TT^uXg (M In 2/2fc B 0) 1/2 , (77) 

which is an inhomogeneous coherence half-length. The small value of s is the main problem of the technique we 
consider here. 

The situation becomes more complicated for the strong acceleration regime, where one has to include the phase 
caused by the non-linear part of the potential [see the second term in the brackets of Eq. (pQ)]. The contribution to 
this phase arising from the (n,m) term of the potential expansion (34t) is given by 
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= (u /u) 2n+1 O , (78a) 
0o = 2 1 - 3 " 7 m" +1 c nm (a, b.) hn,2 m e n 5 2m+n e-\ (78b) 

Leading terms here are those associated with (n,m) given by Eq. (f53|), for which we denote phases as (f>io and 020- 
Including only these phases, one finds that the total phase is given by 

<\>tot = (0 2 o + 0ao + c o) (uo/u) + 0io {uq/u) + 020 (u /u) 2n2+1 . (79) 



Since new terms are not proportional to u , one can- 
not choose a compensating phase cO to reduce tot only 
to the desirable phase (72), but one can choose cO to 
offset the most dangerous contribution to the aberration, 
that linear in a. Cancellation of this term occurs when 



0cO = -0C/O - 50io - (2«2 + 1) 02 



(80) 



For this choice one recovers expression (|75|) for the grat- 
ing profile, in which one has to insert the phase shift 
—4010 — 2n202o, and change the parameter ( to the value 
C = (k g z + 1O0io + n 2 (2n 2 + 1) 2O ) o\. 

We next consider the role of the atomic beam's angular 
divergence. If the angle between the initial momentum 
and the (x, y) plane is non-zero, then the atom enters 
the acceleration zone at a non-zero momentum projection 
along z, p in — M u8 , where Q is of the order of the angu- 
lar divergence. For the weak scattering regime at pi n =/= 0, 
one has to shift momentum change in the definition of 
phase (24b) as dp (t) — > p in + 8p (t) . Phases quadratic in 
Pin are the same for both arms of interferometer, while 
the pi n -independent part is analyzed above. Therefore, 
we can consider only the parts linear in p;„, <pDi, f° r 
which, using Eqs. (p4|), one finds 0Di = PinZi (t) /%, 
where (t) is the atomic z-coordinate along arm i for 
Pin = (i = 1 or 2, see Fig. |^b). Evidently 0£>; is a 
Doppler phase. The Doppler phase difference, 



4> D = 2TTK0 b (X e — x)/\ g , 



(81) 



vanishes at the echo point, x e « L (1 + ////), where L is 
the distance between quadrupoles along x axis, // <C / 
is a force in quadrupole /. A cancellation of the Doppler 
phase at the interference plane is a common property of 
an atom interferometer |l7|| . For the quadrupole beam- 
splitter, we prove that cancellation occurs for a finite 
interaction time, while for optical beam splitters, involv- 
ing couterpropagating waves, Eq. (^Tj) is valid only in 
the Raman-Nath approximation. Owing to the Doppler 
phase cancellation, the only requirement for the weak ac- 
celeration regime is that the angular divergence be less 
than the scattering angle, 



h< (////)• 



(82) 



The situation changes for the strong acceler- 
ation regime, again owing to the phase ( J60| ) 
sensitivity to the non-linear part of potential. 



For pin 0, in the quadrupole 77, 8z (t) = 

M~ l (p in (t + r) - f(t + t) 2 /2j . Assuming that the 

change of the atomic position is small, one finds that the 
additional Doppler phase, which is linear in pi n and as- 
sociated with the (n, m) term in the potential expansion 

(fit]), is s iven b y 



inn. — ^ 



(83a) 

K- n (nnCnrnhn-l.lmr 1 S gS' 2 ^ Q 1 ^ . (83b) 



Requir ing t his phase to be small for leading terms in the 
series ( 34b ), one obtains a condition for the angular di- 
vergence, 



< min {(////) i 



2>> 



(84) 



where 6>i = 6 2 ,mi,0 2 = 6>„ 2 ,o- 

One can use expression (fl|) to estimate the thickness 
Ax of the layer along the x-axis where interference oc- 
curs. Requiring 0c < 1, one finds 



Ax < \ g /2Trn9 b 



(85) 



We can also estimate the width of the layers Sx, in 
which one excites atoms to start and to stop the accel- 
eration. When this width is non-zero, the time of ac- 
celeration is not fixed, and the momentum change Ap is 
spread across a range of width f5x/u. This width should 
be smaller than h/s, i.e. 



5x < 8x = 2d\ g / kits. 



(86) 



Knowing Sx one can estimate parameter (|3j) as tozTi = 
<fiuoSx/2d. 

As an example for the scheme described in 
this paper, consider a beam of 87 Rb atoms 
(L = 0, 5" =1/2, 7 = 3/2), initially pumped into the 
state \G = l,m = 0) and having a longitudinal temper- 
ature O = 1/iK. As we explained in the Introduction, 
acceleration occurs near the centers of the quadrupoles. 
Quadrupole / just splits an atom trajectory into two arms 
(see Fig. ^|b). The main part of the atomic momentum 
change Ap is acquired in quadrupoles // and III. We 
present results of calculations for the last quadrupoles, 
where one can expect the most severe restrictions for 
the system parameters. Inside the quadrupoles, the 
atom is accelerated in the states \G = 2, m = ±1) with a 
magnetic moment /i = ±/xb/2. 
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It is not evident in advance what role the different 
types of nonlinear corrections to the atomic wave func- 
tion play. This role depends on the acceleration regime, 
bias field strength, and the quadrupole geometry. One 
notices, nevertheless, that all parameters of the system 
depend on three variables. Instead of e±, £2, £q, one 
can choose any other three linearly independent parame- 
ters. It is reasonable to choose variables which are most 
severely restricted, and consider how large they can be 



such that nonlinear corrections are still small. 

For a weak acceleration regime we choose the coher- 
ence length s, the layers' thickness Sx and the ratio r\ of 
the atom displacement and the beam radius as indepen- 
dent variables, given by Eqs. ( [43| , [77| , poj ), respectively. 
We found that for 77 = 0.1 and quadrupole size a — 1 cm, 
nonlinear corrections do not rise above 10% if s = 40/ito 
and Sx = 15/im. The beam and fields parameters corre- 
sponding to this choice are given in Table |. 



TABLE I. Parameters of the beam of 87 Rb atoms and the magnetic quadrupole that one can choose to obtain a grating of 
A 9 = 100 nm period (case 1) and lOnm (cases 2 and 3): 6 is the half- width of the incident wave packet; d is the half-thickness 
of the region in which acceleration occurs; u is the beam velocity; B' is magnetic field gradient; b„ and B s are the relative and 
absolute strengths of the bias magnetic field; a is the aspect ratio of the quadrupole size; s is the coherence half-length; £1, £2, 
and e q are relative weights of corrections to the atomic wave function; 9 is the scattering angle; 9i and 62 are upper bounds 
for the beam angular divergence arising in the strong acceleration regime; (j>uo, <t>io, and <j>20 are atom grating phases caused 
by the homogenous and nonlinear parts of the potential; Sx is the thickness of the region in which one starts and stops atomic 
acceleration; uozTi is the parameter (^); Pq and Pq are the geometric averages of traveling wave powers one should apply to 
split atoms between m = ±1 Zeeman sublevels and to produce a 7r-pulse on the transition \G = 1, m = 0) — ► \G = 2, m = ±1) , 
respectively [Po and Pq are evaluated using Eqs. (A9, A14) where we put Sz — 0.1 cm and = 2tt x 1GHz]; 77 is the ratio 

of the atom displacement during acceleration and the beam radius, 771 is the ratio of the time of acceleration and the time of 
wave packet spreading, 772 is the ratio of the momentum distribution width and the momentum change during acceleration; 
parameters 773, 774, 775 verify the validity of the slowly varying amplitude approximation (see Appendix Fa). The parameters 0%, 
62, 0io, 020, ?73, 774 are not relevant in the weak acceleration regime. Parameters marked with stars in the Table are chosen as 
independent. All other parameters depend on these. 



case # 


1 


2 


3 


regime 


weak acceleration 


strong acceleration 


strong acceleration 


X g [nm] 


100 


10 


10 


b[fim] 


104 


3.5* 


A* 


d [cm] 


0.94 


0.55 


0.88 


u [m/s] 


21 


18 


21 


B' [Gs/cm] 


79 


1190 


847 


b s 


0.1 


1 


0.54 


B s [Gs] 


116 


1010 


705 


a 


4.8 


1.5 


1.99 


s [fim] 


40* 


3.5* 


4* 


£1 


0.087 


0.088 


0.11 


£2 


0.027 


0.0073 


0.017 


£ q 


0.0098 


0.0022 


0.0016 


e 


2.2xl0" 3 


0.025 


0.021 


max{6>i, 62} 




0.015 


0.0045 


4>uo 


9.2xl0 5 


5.3xlO y 


5.2xl0 B 


4>w 




6.8 


11 


020 




1.4 


4.0 


Sx [/jm] 


15* 


10* 


14* 




733 


4800 


4200 


P Q [mW] 


0.15 


0.13 


0.15 


[mW] 


0.12 


0.11 


0.12 


V 


0.1* 


20 


24 


Vi 


1.9xl0~ 4 


0.11 


0.12 


V2 


3.1xl0" 4 


9.1xl0" 4 


8.0xl0" 4 


V3 




0.050 


0.43 


Vi 




0.0075 


0.018 


V5 


0.18 


12 


37 
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In the strong acceleration regime it makes sense to 
choose a grating target period X g = 10 nm. As inde- 
pendent variables, we choose the coherence half-length s, 
the beam radius b, and the layer thickness Sx. To achieve 
large values of the parameter 77 and yet small weights for 
the corrections we choose s — b = 3.5/oti and Sx = 10 fim. 
Data for this case are given in the second column of Table 
H A slightly better situation arises if one chooses condi- 
tions, in which nonlinearities quadratic in space, second 
order and fourth order in time vanish. These conditions 
arise for a w 1.99 and b s w 0.54, which are roots of Eqs. 
(|7|, ||). Data for this case are given in the third column 
of Table @. 

Beam splitters can operate also in a pulsed regime. If 
the whole layer 

M < d (87) 

is illuminated by a short raman pulse, atoms are split 
between to ^ Zeeman sublevels and start to acceler- 
ate. A second, time-delayed pulse stops the acceleration, 
and produces two groups of states in the to = sublevel 
with different momenta. When these groups recombine 
at the interference plane, a pulsed atom grating is gen- 
erated. The pulsed grating can be repeated with some 
repetition rate. The pulse regime will have restricted 
application in lithography, because in the time intervals 
between the pulses the flow of atoms will continue pro- 
ducing a uniform background. However, this flow could 
be blocked, by placing a beam stop between the two in- 
tcrferometric arms. The great advantage of the pulsed 
regime is that the time of acceleration is the same for 
all atoms and, therefore, the grating period X g becomes 



velocity-independent. It allows one to relax the severe 
requirements for longitudinal cooling. We will consider 
the pulsed regime in more detail in a future publication. 

In this article we have analyzed the role of the longitu- 
dinal degrees of freedom using a perturbation theory in 
the scattering angle 9. To our knowledge only two other 
articles address this problem jl9],[20) for a finite value of 9 
in an atom interferometer consisting of a set of spatially 
separated, resonant traveling waves. The consideration 
in [p0| that assumes the edges of the field envelopes are 
shorter than an atomic de-Broglie wave length, while a 
quasiclassical approach was used in |l9") . Notably, we 
failed to solve the relevant Schrodinger equation for scat- 
tering from a large angle beam splitter. However, we 
can stress that, when seeking sinusoidal atom gratings 
of a period smaller than the optical wavelength, but still 
larger than the de-Broglie wavelength, our perturbation 
theory is sufficient. If 9 ~ 1 (AdB ~ A) the atom wave 
packet will be broadened in momentum space, which 
would destroy the sinusoidal shape of the atom grating 
and defeat the purpose of our method. 
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APPENDIX A: GROUND STATE DRIVING BY LINA.LIN POLARIZED FIELDS. 



Consider an atom's interaction with a field consisting of two resonant waves propagating along the y-axis, 

E (r,t) = C (x, z) 2_, \^i e ] ex P ( — + iky) + c.c, 

J'=l,2 



(Al) 



where Ej, ej, Qj, k, £ (x, z) are the field amplitude, polarization vector, frequency, wave vector, and envelope function. 
The envelope function £ (x, z) is the same for both fields. For linUin polarized fields one can choose ei = z, e2 = x. 
We assume that Q (x, z) is centered along the atomic beam trajectory, has a small width along the x-axis {Sx <JC d) , 
and a large width along z-axis (8z ^S> b) . The length scales b and d are the atomic beam radius and the acceleration 
zone length, as, respectively, defined in the our paper. 

When field detunings from resonance are larger than the excited state decay rate, the atomic ground state amplitudes 
"0Gm (G and to are the total moment and magnetic quantum number) evolve in the atomic rest frame (x = ut) as |Tst | 



iipG'm' = (G'm \V\ Cm) tpGm, 



(A2) 



where 



(G'm'\V\Gm) 

exp 



jj' 



L0 G'm'Gm i 



A 



(#') 

G'G 



G' K 
to' Vk 



G 

— TO 



Cm") 



(A3a) 
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rC"') 



G'm' ,Gm, 



— Slj — \lj> — tUG'm'; 



=(-i> 



/f+G' + J G / + ,/ G + ./ff +/ 



(A3b) 



(i) ( (*') 

■^JhJg \ a JhJg 



A 



(i) 



V(2A+1)(2G' + 1)(2G+1) 



1 1 



CO?') = (-1)" eie{_;,V2K+l 



J G , I G 
G K J a 

1 1 K 

v v' —Vk 



(A3c) 
(A3d) 



UG'm',Gm IS the G 



Gm transition frequency, Jq and Jh are the electronic angular momenta of the ground and 



excited state manifolds, / is the nuclear spin, X^P H j G — (Jh \ \d\\ Jq) EjC {ut, 0) /2h is the Rabi frequency associated 
with field j, e\, is a spherical component of the polarization vector e^, (...) and {...} are 3J- and 6J-symbols, and we 



have assumed that the detuning Aj G is larger than the Zeeman and hyperfme splitting of the excited states. To be 
specific, we put Jq = 1/2 and Jh = I = 3/2, which corresponds to the D 2 line in 87 Rb. 

Consider first the case where the fields are tuned to the two-photon transitions between Zeeman sublevels of the 
G = 1 manifold, i.e. O2 = Q\ + ujz, where loz = Wn-io = <^io-i — i- The wave function evolves as 



#ii=Xo^n + (x-X*e M ) V10, 

#10 = (x* ~ Xe- 2luJzt ) V11 + (X - X*e 2 ^ zt ) ^1,-1 + Xo^io, 
#1-1 = (x* ~ x e ~ 2lu)zt ) "010 + Xoi>\ 



(A4a) 
(A4b) 
(A4c) 



where xo 



-3- 1 



A { ^ ] (0) + A {22) (0)j and x = 2- 3 / 2 3~ 1 / 2 Af 1 1) (1) are the ac-Stark shift and effective Rabi fre- 

(A5) 



quency associated with transitions between Zeeman sublevels. For a large Zeeman splitting 

^ZTi > 1, 



one neglects the rapidly oscillating terms in Eqs. (A4) to find the wave functions after the field pulses: 

ip^ = i exp [-iA ] {[1 + cos (A)] tfj n - [1 - cos (A)] - n/2sin (A) V^o} > ( A6a ) 

ipf = exp [-«A ] i2 1/2 sin (A) + ipZi) + cos Aip^A , (A6b) 

^£-1 = - exp [-iA ] - (1 - cos(A))7/) n + (1 + cos(A)) - «\/2sin(A) Vfo > ( A6c ) 

where V'g™ on the right hand side are the initial values of the atomic wave function amplitudes, Ao = dtiXOi 

A = 2 1 / 2 J_ dtix are field areas, and x is assumed to be real. If {ipxi, V'ioi V'i'-i.} = {0> lj 0}j ° ne splits 100% of the 
atoms between Zeeman sublevels m = ±1 using a ^-pulse, 



2 



For fields of Gaussian profile, 

C O, z) = exp (-2 (x/Sx) 2 ~ 2 (z/Sz) 2 ^ , 

1/2 

one finds that the geometric average Pq = {P1P2) of the field powers Pi and P2 is given by 



Po 



(2tt) 3/2 hm e c 2 /e 2 Xj H j G f ( J G , J H ) 



(2) 



(A7) 



(A8) 



(A9) 



where m e is the electron mass, \j h j g and / (Jd Jh) are the wavelength and oscillator strength associated with the 
excited-ground state transition. 

Atoms in |m = ±1) sublevels start to accelerate in an inhomogeneous magnetic field. To stop this acceleration at a 
later time, one needs to return the atoms back to the |m = 0) state. Inserting initial conditions ip^ x = 1 (or V'iT-i = 1) 
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in Eqs. (|A6j) one sees that one can return at most half of the atoms to the \m = 0) state, again using a ^-pulse. The 
other half remains split between the |m = ±1) states. 

This loss of atoms can not be avoided in the quadrupole / (see Fig. ||), if one operates in the weak acceleration 
regime. However, for the strong acceleration regime or for quadrupoles // and one can use different fields along 
different arms of the interferometer and employ another hyperfine sublevel to achieve a 100% exchange between 
accelerated and non-accelerated Zeeman sublevels. 

For example, to start the acceleration in quadrupole II, one chooses the field frequency difference 



^2 — Ol = U>21, 



(A10) 



such that under condition (A5), only the transition between sublevels \G — l,m = 0) and \G = 2,m = 1) occurs. The 
wave function amplitudes of this two-level system evolve as 



#21 
#10 



xiV^i + xV'io, 

X^2l +Xo^io, 



(Alia) 
(Allb) 



< 3 (0) 



where X i = ^ 1/2 (A? (0) + A% 2) (0) J , Xo 

all atoms between the sublevels, one needs the ac-Stark shifts to be equal, Xi 
fields' powers has to be chosen as 



+ A {22) (0)) , X = 2- 3 /25-i/2^(2i) (!) _ To transfer 
Xoj which means that the ratio of the 



P2/P1 = 

This ratio is positive only for negative detuning, 



-1 - 2w 21)10 /a£ )1 



(A12) 



-2cj 2 i 



< A 



(i) 



< 0. 



To obtain a 100% transfer between the levels, one should apply a 7r-pulse, for which dt \ X \ = n/2. Thi s condition 
is an equation for the geometric average of the field powers. Combining this equation with the Eq. ( |Al2j ), one finds 
the powers, 



P2 = Pn 



2w 2 i,io/Aj2 i1 
2cj 2 i,io/Aj f | jl 



-1/2 
1/2 



where 



P' - 
^0 — 



A^I 2 hm e c 2 /i 1 / 2 e 2 Xj H j G f ( J G , J H ) 



uSzA 



(2) 



(A13a) 
(A13b) 

(A14) 



APPENDIX B: JUSTIFICATION OF THE SLOWLY VARYING AMPLITUDE APPROXIMATION. 



In this Appendix we determine conditions under which it is valid to assume that the amplitude ip (z, t) in Eq. 
( |i"4| ) varies slowly, i.e. the operator Q given by Eqs. (|2^) leads to small corrections to the zero-order approximation 
solution (|3"(i|). To find these conditions, one has to include the time dependence of the force and the homogeneous 
part of the potential. At small times, \t\ <C r a , these are given by 



/(*) «/(i + £* 2 A- a 2 ), 

U (t) =U (1 + ut 2 / T l) 



(Bla) 
(Bib) 



where v ~ 1 and £ is given by 



i = 2 [3&^ (a 2 - 1) (a 2 + l) 2 - 16a 2 &7 2 (a 2 



(B2) 



for a potential produced by a magnetic quadrupole. We found that if (i) one can neglect the first term in Eq. (26a) 
and (ii) approximation (^) is valid, then the slowly varying amplitude approximation is valid if the correction p£) 



is of a small relative weight 



e q < 1) 



We now prove assumptions (i) and (ii). 
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We start from the first term in Eq. (26a). During the acceleration it produces corrections to the wave function of 
relative weight J^ r dtj^ — f* d ~ ^ and Auhq/Mu 3 , where Au is a typical change of the velocity (|l3|). Under 
condition (|||), using the estimates 4f ~ (U /Mu 2 ) S 2 , U ~ fa, 

q « /9 P0 ~ (B3) 

421) , one finds that the weights are small (~ 68 <C 1 and 9 2 [3 <C l) , and one can eliminate the first term 



and Eqs. 
in Eq. (|6e 



Now, consider the remaining part of Eq. ( 26a ). For the operator q one finds 

q = d t - fd Po -i(<t>- f(t) 8z/h + p 8z/h) . 



(B4) 



Calculating <p and 8z using Eqs. (E4, Bl) one finds 



, u u? {(9t ~ fd Po f + i (Mh)- 1 [fSp - p f - e/ 2 
■ 2 (poSp + e/ 2 r 4 r- 2 ff ) (d t - fd Po )] ~ (Mhy 2 (p 8p + 2C/Vr; 2 gf) 



(B5) 



where g= -3 (t/r) 4 - 8 (t/rf - 6 (t/r) 2 + 1 /12, .g' = rdg/dt. 

Since we have used the operator Q above only for the evaluation of the corrections, it is sufficient to consider 
the operator Q acting only on the unperturbed wave function (PQi~t)- For weak acceleration, when Wo (po>i) — 
{poi ~ T ) j one can eliminate the time-derivative in Eq. (B5). For strong acceleration, owing to the phase associated 
with the second term in brackets in Eq. (|60|), ^o(po>i) = exp — £ J_ dt\Ux [Sz (ti) ,ti] Wo (.Po,~ T ) , the time- 
derivative leads to the factor 



d t ~ Ui [Sz (t) , t] /%, 



(B6) 



The leading independent terms in the U\ expansion (34b) are associated with the (2, mi) and (ri2,0) elements in 
Table (|39j) . Retaining only these terms and assuming that 8z (t) ~ Az and t ~ r, one obtains the estimate, 

d t ~ /Az 2 (ftar 1 (5 2mi + /Az^rV 1 -" 2 . (B7) 



Comparing this result with estimate (B3), one finds that the time-derivative can still be eliminated even for strong 
acceleration if parameters 



V3 



m = (soy 12 /r 1 



■i<\) 

(B8b) 



are small 

ing for e 

factors r}i, 171772, <5 2 ?7i, rj 5 , vh ViVs, where 



Assuming for estimates that 8p ~ /r, one finds that the weights of other contributions to q 2 differ from (fd Po ) by 



775 = os 3 ?- 1 . 



(B9) 



From Table Q in Sec. IV, one sees that all these factors are smaller than unity, except the factor 775, which is large for 
the strong acceleration regime. The contribution of the or der of 775 is proportional to t he p arameter £. Since parameters 
f and C21 vanish simultaneously [compare Eqs. ( |38a| , B2)], one can exclude in Eq. ( |B5| ) all terms containing £ if the 
aspect ratio of the quadrupole a and the relative bias field strength b s are chosen to satisfy Eq. (|57j). Therefore, for 
all cases under consideration, the expression ( [l8|) provides the main contribution to the correction associated with a 
slowly varying amplitude approximation. 
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